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The effective ^-matrix from conductance data 
in a quantum wave guide: Distinguishing the indistinguishable 
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We consider two different stationary random processes whose probability distributions are very 
close and indistinguishable by standard tests for large but limited statistics. Yet we demonstrate 
that these processes can be reliably distinguished. The method is applied to analyze conductance 
fluctuations in coherent electron transport through nanostructures. 
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The quantum-mechanical scattering by complex sys- 
tems has been a problem of long-standing interest. This 
problem is of a great importance in nuclear, molecular 
and mesoscopic physics [J, B 0, ■ The scattering pro- 
cess can be often analyzed in terms of two distinct time 
scales: (i) a prompt response due to direct processes 
and (ii) a time-delayed scattering mechanism. In nu- 
clear physics the time-delayed mechanism is associated 
with the formation of an equilibrated compound nucleus 
0, 0| • Similarly, time-delayed processes have been en- 
countered in the study of coherent electron transport 
through nanostructures 0, 0, 0, El ■ The two distinct 
time scales introduce two different energy scales. Namely, 
the scattering amplitude of direct processes is a smooth 
function of energy and it can be taken as an energy aver- 
aged S-matrix. The scattering amplitude of time-delayed 
process, on the other hand, is an energy dependent func- 
tion which fluctuates rapidly around zero. Knowledge of 
the energy averaged ^-matrix allows one to construct a 
probability distribution (PD) of the full S'-matrix 0,0,0. 

S'-matrix for ballistic electron scattering off nanostruc- 
tures can be calculated numerically if the geometry of 
the nanodevice, as well as the shape, size and number of 
propagating modes of the leads, and spatial distribution 
of disorder are known precisely Then, having calculated 
the energy dependence of the ^-matrix one can directly 
separate it into an energy averaged direct reaction com- 
ponent and a fluctuating one. In experiments, the de- 
tailed information about the conducting device may not 
be always available. For example, wall imperfections and 
unknown spacial distribution of disorder does not allow 
to evaluate the contribution of direct processes. Disor- 
der can also block leads in an unknown manner reduc- 
ing the number of propagating modes. This motivates 
us to consider the "inverse problem" : What information 
about relative contributions of direct and time-delayed 
processes and number of propagating modes in the leads 



can be obtained from the analysis of the transmission 
energy fluctuations provided the above mentioned char- 
acteristics of the nanodevice are unknown? 

For concreteness we consider the transmission between 
horizontally oriented leads for the two types of geometry 
of wave guides (WG) (Fig. 1). For experimental pur- 
poses, we refer the reader to reference |9| for a possible 
implementation of this system. The left leads of both 
WG can accommodate a single propagating channel only. 
For the WG (b) the right lead also accommodates only 
one channel (N = 1), while for WG (a) it has N > 2 
open channels. In addition, both WG (a) and (b) have 
N a ^> 1 and Nf, ^> 1 open channels shown by arrows 
pointing down at the bottom of WG with N a ^> N. Dis- 
order is modeled by scatterers (shown in Fig. 1 in white) 
having infinitely high potential walls. 

Transmission between the horizontally oriented leads 
for the (a) and (b) WG in Fig. 1 is given by 



A' 



T(E) = J2\Si j (E)f 



(1) 



Here Sij (E) are the S- matrix elements for scattering be- 
tween left leads "1" and right leads j, and E is electron 
energy. For the WG (b) in Fig. 1 the sum (1) con- 
tain only a single term (N — 1). We use the decom- 
position Sij(E) =< S\j > +5Sij(E), where < Sij > 
are energy averaged ^-matrix elements, and 5S%j(E) are 
energy fluctuating ones. We chose the spacial distribu- 
tion of disorder such that (i) JY | < Sij > | 2 <C 1, and 
(ii) transmission between the horizontal leads and open 
channels at the bottom of both (a) and (b) WG in Fig. 1 
is mostly due to time-delayed processes. Therefore, the 
overall amount of direct processes for the transmission 
between different channels is much less than that for the 
overall amount of time-delayed processes. Also, due to 
the presence of disorder, the dynamics of the classical 
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FIG. 1: The electron probability distribution inside the wave 
guides for an arbitrary chosen electron energy. The arrows 
show incident unit current from left and current leaving the 
system from the right. The white circular regions represent 
infinite potential barriers. The wave guide (a) does not show 
direct processes and (b) has direct processes (see text). 

counterpart in the interaction region of WG (a) and (b) 
in Fig. 1 is chaotic. Since the total numbers of open 
channels, A" + N a + 1 and N + 2, are much greater than 
unity, we are in the regime of Ericson fluctuations when 
electron resonance states within the WG are overlapping. 
Then each individual fluctuating S'-matrix element can 
be considered as a Gaussian random process. Namely, 
real and imaginary parts of each individual SSij(E) in 
sum (1) are distributed by a Gaussian law, are uncorre- 
cted and have the same dispersions 0, |{| . 

The PD for y = T(E)/ < T(E) > (so that < y >= 1) 
is given by 0, 

P N (y) = Nil-yar^y/y^-V' 2 
exp[-N(y d + y)/(l-y d )} 
7 A r_ 1 [27V(y £i2/ ) 1 /2/(l_y £i )]. (2) 

Here I/v_i(A) is the modified Bessel function of (N — 1) 
order and y d — J2 3 ; I < Slj > IV < T(E) > is tne rela- 
tive contribution of direct processes to the total trans- 
mission. The PD of Eq. (2) is exact if all 5Sij(E) 
with different j are uncorrelated and have the same dis- 
persions. Otherwise this PD is an approximation for 
which N should be understood as an effective number 
of independent channels N e ff < N. It is given by 
the normalized variance of ST(E) = \3Sij(E)\ 2 '- 
< ST(E) 2 > j < 5T{E) > 2 -1 = 1/N eff . 

For the WG (b) in Fig. 1 one has to put N = 1 in 
Eq. (2). For y d — 0, Eq. (2) takes the form of a \ 2 - 
distribution with 2N degrees of freedom 

P N (y) = NFiN)]- 1 ^)"- 1 exp(-Ny), (3) 




FIG. 2: The probability distributions of transmission. The 
solid lines, with decreasing thicknesses, show the theoretical 
curves for N = 1, N = 5 (Eq. (2)) and N = 10 (Eq. (3)) 
channel cases respectively. The histogram is constructed from 
the numerical data for TV = 1 and yd = (0.9) 1 ' 2 , as an exam- 
ple. The inset shows a blow up of the transmission probability 
distributions demonstrating that the difference between them 
around their maxima does not exceed 4%. 



where V{N) is Gamma function. Note that the variance 
of y is given by < y 2 > -1 = (1 - y d )/N @- 

As an example, we consider two possible cases for WG 
(a) in Fig. 1: (i) N = 10 and y d = 0, the latter due to 
complete blocking of the direct paths between the leads 
by disorder, and (ii) N = 5 and y d — (0.8) 1 / 2 . Also sup- 
pose that, for WG (b), for which N — 1, y d — (0.9) 1 / 2 . 
For the above three sets of A^ and y d values the variance 
of conductance y is the same and equals 0.1. The ques- 
tion is: Can we distinguish these three cases by analyzing 
PD of y in Eq. (2)? 

In Fig. |2 we plot the PD for the three sets of N and 
y d . The three distributions are very close. We performed 
a x 2 test and found that, even for as many as 28000 in- 
dependent realizations of y for each of the three cases, 
the distributions are indistinguishable at a 99% confi- 
dence level. And, to the best of our knowledge, there is 
currently no any other statistical test which would allow 
to distinguish between the three different stochastic pro- 
cesses. Yet, in what follows, we show that the problem 
of distinguishing between these random processes can be 
solved. 

We will refer to the j/ as j/i, where the index i stands 
for different independent realizations of the process. Let 
us transform to new random variables = (y; + yj)/2, 
r ij = (Vi ^ Uj)/^ s ij with i ^ j. One can now easily find 
the joint probability distribution K(s, r). The analytical 
calculations are formally similar to those of Ref. [l(| for 
the analysis of the correlation between cross section and 
analyzing power in the regime of Ericson fluctuations in 



3 



oT 0.05 
v 0.048 





N ,h =1 ' 




OT^I 




- - N ,h =5 




— tr»=5 . 




... N * =10 




^""=10 ■ 




T 


1 ■* ™ 



s 

FIG. 3: The plot of < r 2 (s) > vs. s corresponding to the 
same TV (shown in inset) and yd values as those used for the 
analysis of the probability distributions in Fig. |U Dotted, 
dashed and dotted-dashed lines are theoretical predictions 
(Eq. (4). 



nuclear reactions. Without presenting the explicit result 
for K(s, r) we point out that, for yd — 0, s and r are 
statistically independent. On the contrary, for yd ^ 0, s 
and r correlate, the bigger yd the stronger this correlation 
is. 

We calculate < r 2 (s) >: 

< r 2 (s) > = [1 - / 2iV+ 3(A)/72JV-i(A)]/(2iV + 1), (4) 

where A = 4N(yds) 1 ^ 2 /(l — yd), and take into account 
the constraint (1 - yj)/N = 0.1. In Fig. 01 we show 
< r 2 (s) > for the three cases. We have generated 28000 
realizations for each case and used all possible i < j 
combinations for s.y and nj . The data in Fig. [3] are 
obtained by dividing the s%j into 5 bins for each case 
and taking the average of r 2 (s) and s for each bin. It 
is evident from Fig. that we can distinguish the three 
cases. The fluctuations around the theoretical curves are 
determined by the number of realizations. Error bars are 
calculated numerically from standard deviations of 100 
data sets each consisting 28000 y realizations. We found 
that it is possible to distinguish the N — 1 and N = 10 
processes in 97% of cases for 210 realizations of y. This 
drops to 85% for 140 realizations. 

We apply the new method to distinguish between (a) 
and (b) in Fig. 1 from the conductance fluctuation 
data. We implemented a finite element solution of the 
Schrodinger equation 0, |n|. Calculations were per- 
formed on the energy range which allows 5 propagating 
channels in the right lead for WG (a) in Fig. 1. Both 
left leads for WG (a) and (b) as well as the right lead for 
WG (b) in Fig. 1 accommodate a single channel. In ad- 
dition, for both (a) and (b) WG in Fig. 1, there are 20-25 
propagating channels from the lower part of the devices. 
A unit current is directed from left. We are interested 
in the transmission probability from the left to the right 
side. Spacial disorder distribution for WG (a) in Fig. 1 
was chosen to block the direct paths between the left and 
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FIG. 4: The plot of < r 2 (s) > vs. s for the wave guides 
shown in Fig. Q The dashed and dotted-dashed lines are 
theoretical estimates (Eq. (4)) and thin and thick solid lines 
are numerical calculations for TV = 1 (with direct processes) 
and for N — 5 (without direct processes) cases respectively. 
The inset shows conductance distributions (lines) for (a) N = 
1, with direct processes (Eq. (2)), and (b) N — 5, without 
direct processes (Eq. (3)). Histograms are constructed from 
the numerical data. 



right leads, i.e. to maximally suppress direct processes. 
This was confirmed by S'-matrix numerical data, which 
yielded negligible absolute value of the energy averaged 
S-matrix, and by the fact that the normalized variance 
of conductance fluctuations is close to 0.2. For WG (b) 
in Fig. 1 the extension of leads into the cavity is adjusted 
to have relative contribution of direct process yd ~ 0.9 so 
that the normalized variance of conductance fluctuations 
is close to 0.2, i.e. to that for WG (a) in Fig. 1. For both 
(a) and (b) WG in Fig. 1 we used different disorder dis- 
tributions to increase the statistics. The overall number 
of independent realizations of transmission values y for 
each type of geometry in Fig. 1 was around 150. 

In Fig. 4 we present PD of y for the (a) and (b) WG 
of Fig. 1. x 2 test does not allow to distinguish the two 
PD. Yet, from Fig. 4 one can see that for the N 
1 case, in the presence of direct processes, < r 2 (s) > 
strongly increases when s decreases. On the contrary, 
for the N = 5 case, when direct processes are absent, < 
r 2 (s) > for s < 1 is close to a constant. We have checked 
that even for each single configuration of disorder, when 
we have only about 30-40 independent realizations of the 
transmission values for each (a) and (b) WG, < r 2 (s) > 
increases always noticeably faster for the N = 1 case with 
direct processes as compared to the N = 5 case without 
direct processes. On the other hand, the rise of < r 2 (s) > 
when s decreases of s for N = 1 with direct processes 
is much stronger than that predicted by Eq. (4). In 
order to find a possible reason for this we analyzed the 
PD of fluctuating S'-matrix elements. We found that (i) 
SSn(E) are not distributed isotropically in the complex 
plane, (ii) the PD of real and imaginary parts ofSSn(E) 
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are not Gaussian, and (iii) real and imaginary parts of 
SSn(E) are correlated. Therefore, the reason for the 
discrepancy between the data and theoretical line in Fig. 
4 for N = 1 is that the conditions to derive Eq. (4) arc 
not met. Yet, the distribution for transmission is not 
sensitive to the deviation of the PD of SSu(E) from its 
isotropic Gaussian distribution in the complex plane and 
correlation between real and imaginary parts of SSu(E). 

In conclusion, we have suggested a method to dis- 
tinguish between different stationary random processes 



whose PD are very close and indistinguishable by stan- 
dard tests for large but limited statistics. The method 
has been applied to analyze conductance fluctuations 
in coherent electron transport through WG. We found 
that the method proposed here, unlike the analysis of 
PD of transmission, is sensitive to (i) the deviation of 
PD of fluctuating transmission amplitudes from isotropic 
Gaussian statistics in the complex plane, and (ii) correla- 
tions between real and imaginary parts of the fluctuating 
transmission amplitudes. 
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